Non-equilibrium Pattern Modes Extracted from Experimental Data 
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We describe a method to extract from experimental data the important dynamical modes in 
spatio-temporal patterns in a system driven out of thermodynamic equilibrium. Using a novel 
optical technique for controlling fluid flow, we create an experimental ensemble of Rayleigh-Benard 
convection patterns with nearby initial conditions close to the onset of secondary instability. An 
analysis of the ensemble evolution reveals the spatial structure of the dominant modes of the system 
as well as the corresponding growth rates. The extracted modes are related to localized versions 
of instabilities found in the ideal unbounded system. The approach may prove useful in describing 
instability in experimental systems as a step toward prediction and control. 

PACS numbers: 47.20.-k, 47.52.+j, 47.27.Cn 



Identification of instabilities plays a crucial role in our 
understanding and description of the the dynamics of 
many nonlinear physical, biological, and chemical sys- 
tems driven out of equilibrium [1 . In particular, quanti- 
tative description is essential for predicting and/or con- 
trolling the evolution of such systems, with weather pre- 
diction being a prime example. While linear stability 
analysis of global disturbances in an idealized, infinite 
system may provide a description of dynamics and pat- 
tern selection, this approach fails for imperfect patterns 
(e.g., far from onset) and in strongly confined systems. 
Despite recent numerical advances [2 -4 in computing the 
spatial structure and dynamics of localized disturbances 
in weakly chaotic patterns, no general approach has been 
developed for extracting such dynamical information di- 
rectly from experimental measurements. 

In this Letter, we present such an approach, illus- 
trating how dynamical degrees of freedom can be ex- 
tracted from experiments conducted on the prototypical 
Rayleigh-Benard convection (RBC) system. Specifically, 
we determine the spatial structure and evolution of the 
dominant dynamical degrees of freedom by analyzing the 
response of the system to an ensemble of localized per- 
turbations about the straight roll state. In each case, 
we find the dynamics are dominated by a small num- 
ber of spatially-localized modes. We also observe slowing 
down of the dynamics and thus quantify the distance to a 
particular instability boundary in terms of perturbation 
lifetimes. The spatial structure of the extracted modes is 
found to be consistent with the classification of secondary 
instabilities of spatially infinite perfect patterns. 

The convection experiments were performed with a 
layer of sulfur hexafluoride (SFe) gas of depth d = 
700 ± 10 jam compressed to 14.5 ±0.1 bar [5 . The gas 
layer was confined laterally to a 25 mm x 15 mm rectan- 
gular region by filter paper sidewalls chosen to match 
the fluid conductivity. The layer was bounded from 
above by a water-cooled sapphire window and from below 
by a carbon-disulfide (CS2)-cooled zinc selenide (ZnSe) 
window; the temperatures of these windows were reg- 
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FIG. 1: Experimental setup includes a shadowgraph visual- 
ization as well as a system for optical actuation. Inset are 
observed responses to local perturbations at different param- 
eter values. These disturbances represent the local version of 
the skew- varicose and cross-roll instabilities, respectively, of 
the unbounded system. 



ulated to ±0.05 °C. Experiments were performed with 
T = 24.00 °C and temperature differences of 2.92 and 
4.80 °C. Convection was visualized using the shadow- 
graph method by illuminating from above and imaging 
the light reflected from the ZnSe surface at the bottom 
of the gas layer. 

SFq is a greenhouse gas; we use this property to opti- 
cally apply controlled thermal disturbances both to ma- 
nipulate the global convective flow as well as impose lo- 
calized perturbations. An IR beam from a CO2 laser 
(10.6 /im wavelength) is focused to approximately 200 
/im in diameter and strikes the gas layer from below (af- 
ter passing through the IR transparent CS2 coolant and 
ZnSe window). The beam is steered by two computer- 
controlled, gold-plated servo mirrors that rotate about 
orthogonal axes, providing the ability to direct the IR 
light toward any point over the cell domain. Fig. [T] shows 
the experimental setup. The extinction length of the 



2 



beam in the SFe is < 10 /im 0, less than 2% of the ceh 
depth, so the absorbed beam induces a highly localized 
heating that takes place very near the bottom of the gas 
layer. Software developed in-house works with a com- 
mercial program (LaserShow2000) to synchronize laser 
power with mirror rotations. The strong absorption and 
rapid scanning results in the ability to manipulate the 
flow on a time scale much faster than the typical dynam- 
ical time scale (vertical thermal diffusion time of 2.7 s). 
This technique improves on previous attempts to manip- 
ulate convection patterns by optical actuation [71 18| . 

A pattern of straight rolls with wavenumber q was es- 
tablished by sending pulses of laser light into the cell at 
desired hot roll locations. To minimize sidewall effects, 
the pattern was imposed in the central portion of the 
cell, 1-2 wavelengths from either of the sidewalls in the 
wavevector direction. This provides room for 7-9 interior 
pairs of straight rolls. 

Secondary instabilities of the straight roll state are pre- 
dicted to define the boundaries of a stability (Busse) bal- 
loon in (e, q) space at fixed Prandtl number; crossing one 
of these instability thresholds results in a re-organization 
of the pattern [9 . Here, e = {R — Re)/ Rc is the reduced 
bifurcation parameter that measures the distance from 
onset of the primary instability of the purely conducting 
state. Fig. |2] shows the Busse balloon for the spatially 
unbounded system at the conditions of our experiments 
(Pr = 0.84). The Busse balloon is conventionally cal- 
culated assuming an infinite system; we may expect the 
thresholds for the localized instabilities in a bounded sys- 
tem to take place at slightly different parameter values. 
Nonetheless, the balloon provides a useful reference as 
different areas in the parameter space are visited, and 
the dynamics of the dominant modes of the system can 
be expected to depend on the distances between a given 
(e, q) point and the various instability boundaries. 

Pattern manipulation allows us to use both e and q 
as control parameters. The initial pattern wavenumber 
was chosen to be within the stable band. Following the 
inital imprinting, actuation is restricted to only the two 
outer- most rolls between which the straight roll pattern is 
confined. The pattern wavenumber is adjusted by moving 
the positions of these outer rolls toward or away from one 
another; the interior pattern equilibrates on a timescale 
of approximately 10 for our domain size. A closed- 
loop feedback algorithm constantly analyzes images for 
deviations of the two boundary roll positions from desired 
locations and adjusts laser power accordingly. 

Sufficiently close to a particular instability boundary, 
the spatial modes corresponding to that instability are 
weakly damped and can therefore be excited by small 
perturbations to the base state (straight rolls). Each 
perturbation takes the form of a brief (about 100 ms) 
well-localized laser pulse directed to a single location in 
the convection pattern. We estimate that the applied 
heating results in an initially axisymmetric disturbance 
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FIG. 2: Stability balloon at Pr = 0.84 showing various insta- 
bilities (CR = Cross Roll, ECK = Eckhaus, SV = Skew Vari- 
cose, OSC = Oscillatory) of the straight roll state along with 
the parameters at which the experimental data were taken for 
the three data sets. Di = (0.60, 2.85), Dn = (1.50, 2.40), and 
Dm — (0.60, 2.20). The inset plot shows the disturbance de- 
cay rate (s~^) as a function of wavenumber q at fixed e = 0.60. 



of diameter d/2. 

We first probed the pattern response to perturba- 
tions for different q near the high- wavenumber instability 
boundary, with fixed e = 0.60. Upon perturbation, the 
two hot rolls appear to bend toward one another and 
then relax back to their initial locations; a snapshot of a 
typical response is shown in the inset of Fig. [l] The dis- 
turbance decay is slower at higher which suggests using 
the perturbation lifetime as a measure of distance from 
instability. The perturbation lifetime is measured most 
easily from the apparent local roll separation as a func- 
tion of time. This signal typically displays a large spike 
immediately after the perturbation, followed by exponen- 
tial decay to the original value; analysis is restricted to 
this period of linear decay. 

We find that the disturbance lifetime increases linearly 
with the pattern wavenumber over an order of magnitude 
change in decay rate (see the inset in Fig.[2|. The decay- 
rate distribution shows the expected slowing down of the 
dynamics near the instability and indicates the critical 
value qc — 3.01 for the localized skew- varicose instability. 
Note that this value is slightly smaller than that of the 
global instability predicted from analysis of an infinite 
domain ((^c = 3.15), which reflects the effect of spatial 
localization. The low amount of scatter in the decay- 
rate plot illustrates the high degree of reproducibility of 
the imposed perturbations. 

In order to excite all dominant localized modes, per- 
turbations were applied at a grid of equally-spaced loca- 
tions across a wavelength of the pattern. All other per- 
turbations are related to this set through the symmetries 
of the system (translational invariance in the direction 
perpendicular to the wavevector, periodicity in the di- 
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rection of the wavevector), as long as the disturbances 
are sufficiently well-localized and not near the physical 
boundaries. Each set contains perturbations at 12 dis- 
tinct locations; in all experiments, the spatial extent of 
the laser is of order 10% of the pattern wavelength, so 
there is no benefit from a finer partition. 

Shadowgraph images capture the full evolution of an 
imposed disturbance and thus contain snapshots of the 
composite structure of the excited modes over time. Dur- 
ing the perturbation decay, there exists a period of time 
over which the dynamics of the disturbance can be de- 
scribed by a linear evolution operator; we seek to rep- 
resent this evolution operator by computing its matrix 
elements in a subspace spanned by a set of slow modes 
extracted from the shadowgraph images. We first sub- 
tract the stationary straight roll pattern from all images 
after the perturbation. The images are then spatially 
windowed and Fourier filtered. A Karhunen-Loeve (KL) 
decomposition of difference images representing pertur- 
bations 2-5 seconds after each initial disturbance provides 
a set of basis modes. All perturbations are then projected 
onto the subspace spanned by these basis modes. Note 
that the typical implementation of the KL decomposition 
uses time- averaging [10], whereas we employ an ensem- 
ble average over different initial conditions. We limit our 
embedding dimension to a small (usually three or four) 
number of modes that capture 90% of the power of the 
KL eigenvalue spectrum. 

Let us denote the disturbance following an initial per- 
turbation b^, expressed in a low-dimensional basis. Af- 
ter some time T, this state has evolved to b"^. Then 
Uh^ = b"^, where U is the evolution operator. Using an 
ensemble of initial conditions we define 

B°=[b?|b°|...|b°,] 

and similarly, 

[bf I I ••• I hi] 

This gives the over-determined (least-squares) problem 
for the evolution operator UB^ = ^ which is solved by 
U = B^{B^) ^, where the reciprocal of B^ is taken to 
refer to the generalized inverse of the non-square matrix. 

During linear decay, each system eigenmode decays 
at a characteristic rate, so we can also write Uei = 
exp(crT)ei, where are the eigenvectors and the eigen- 
values Xi are related to the growth rates by = 
exp(criT). 

The discrete translational symmetry of the straight roll 
pattern implies that a superposition of a well-localized 
eigenmode with a copy of itself translated over an inte- 
ger number of wavelengths is an eigenmode with the same 
growth rate. Hence, there exist many possible represen- 
tations of each eigenmode. Reduction of the dynamics 
using the symmetries of the base pattern is useful for 
both unbounded and bounded systems. 




FIG. 3: (a) The fundamental dominant and (b) sub-dominant 
mode extracted at Dj. Dashed lines mark the approximate 
locations of the hot rolls of the underlying base state; cold 
rolls lie directly between the hot rolls. 




FIG. 4: (a) The fundamental dominant and (b) sub-dominant 
mode extracted at Dm. 

In particular, all disturbances can be decomposed in 
terms of symmetries related to a pair of hot and cold 
rolls. By defining two symmetry planes, one at the center 
of a hot roll, the other at the center of an adjoining cold 
roll, we can define four symmetric versions of every ini- 
tial disturbance, each even/odd about the hot/cold roll. 
Each of the corresponding four subspaces is invariant: 
disturbances retain their symmetry as they evolve. 

The entire collection of initial and final conditions ex- 
tracted from experiment was decomposed using these 
symmetries, producing four independent ensembles. All 
eigenmodes extracted from the four ensembles are eigen- 
modes of the system, but in the cases when multiple 
eigenmodes share an eigenvalue (growth rate) we elim- 
inate redundant representations by computing the most 
spatially-localized eigenmode structure. 

We estimate the uncertainty in the growth rates to be 
less than 10%. This allows one to group all extracted 
modes and define the fundamental mode as the most- 
localized structure with a particular growth rate. The 
fundamental modes were computed by minimizing the p- 
norm {p < 2) among linear combinations of all modes 
(and their translated copies) in each group. The symme- 
tries imply that all extracted modes in each group can be 
represented as linear superpositions of the fundamental 
mode along with its translated and/or reflected copies. 
We verified that this is indeed the case as such represen- 
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tat ion was accurate, with mutual projection > 0.94. 

The first ensemble was produced at the point Dj of the 
parameter space (see Fig. [2|. The two dominant modes 
are shown in Fig. |3] As the lifetime measurements indi- 
cate, the least-stable mode (ai = -0.13 s~^) tended to 
be excited from perturbations between two hot rolls (at 
the location of a cold roll). Note, however, that while the 
structures excited from these disturbances are even about 
the cold roll, the most-localized representation of this 
mode does not obey that symmetry. The sub-dominant 
mode ((72 = -0.70 s~^) tended to be excited from pertur- 
bations directly at the location of a hot roll. 

A second set of perturbations at high q was performed 
at Dji^ with e increased relative to Dj. Again, two modes 
are extracted; ai = -0.15 and (72 = -0.55 s~^. Mutual 
projections of the modes extracted at Dj and Du indi- 
cate that the spatial structure of the two dominant modes 
remains unchanged (after scaling by the wavelength) be- 
tween these two locations. There is also, in both cases, 
a large separation between the two growth rates. This 
suggests we can identify the dominant mode (Fig. |3^) as 
the one representing the secondary instability at the high 
wavenumber boundary. Its structure is consistent with a 
skew-varicose type instability. 

We also created an ensemble of perturbations to a low- 
wavenumber pattern, at Djjj. As in the other two cases, 
two modes were extracted, with growth rates ai = -0.20 
and (72 = -0.27 s~^. As seen in Fig. [ij the spatial 
structure of the dominant mode does not resemble any of 
the previously extracted modes, while the sub-dominant 
mode resembles the sub-dominant mode extracted at 
both Dj and Du. The closeness of the growth rates 
is consistent with the existence of two low-wavenumber 
instability types of the unbounded system that occur at 
very nearly the same parameter values, namely, the Eck- 
haus and cross-roll instabilities (see the Busse balloon 
in Fig. [2|. The dominant mode is again excited from 
perturbations to a cold roll; this is not surprising, as 
heating of cold rolls tends to reduce the amplitude of the 
saturated state. We find experimentally that sufficiently 
strong perturbations of this kind result in the growth 
of rolls perpendicular to the base pattern. We therefore 
identify the dominant low-wavenumber mode with the lo- 
calized cross-roll instability and the sub-dominant mode 
with the localized Eckhaus instability. 

Further experiments are needed to explore the applica- 
bility of this approach to states exhibiting more complex 
dynamics. One such state, occurring in gas convection 
experiments with Pr ~ 1, is the spatiotemporally chaotic 
state known as Spiral Defect Chaos (SDC) [11 . While 
there exists a bistability between stationary straight rolls 
and SDC [12 over the parameter range of the experi- 
ments reported here, localized instability in the straight 
roll pattern introduces defects which tend to lead to a 
disordered pattern, thus providing a mechanism for the 



transition to chaotic behavior. Additionally, it was deter- 
mined in a numerical study [2^ that the chaotic dynamics 
of SDC are largely driven by the creation/annihilation of 
defects occuring in straight roll regions of the pattern. 
We expect, therefore, that spatially- localized modes are 
dynamically-important in both the transition to and the 
driving of chaotic behavior, suggesting a natural exten- 
sion of our experimental approach to investigations of 
more complex convection patterns. Moreoever, the out- 
lined procedure is general enough to be used in a va- 
riety of other dynamical systems, so long as an appro- 
priate means of system actuation can be developed. In 
addition to being of fundamental interest and of use in 
increasing predictive power, knowledge of the modes of 
instability could be particularly advantageous in system 
control, where small, controlled perturbations could be 
used to guide system dynamics [13 . 
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